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Abstract 

A cryogenic test article, the Generic Research 
Cryogenic Tank, was designed to qualitatively sim- 
ulate the thermal response of transatmospheric vehi- 
cle fuel tanks exposed to the environment of hyper- 
sonic flight. One-dimensional and two-dimensional 
finite-difference thermal models were developed to 
simulate the thermal response and assist in the design 
of the Generic Research Cryogenic Tank. The one- 
dimensional thermal analysis determined the required 
insulation thickness to meet the thermal design crite- 
ria and located the purge jacket to eliminate the lique- 
faction of air. The two-dimensional thermal analysis 
predicted the temperature gradients developed within 
the pressure- vessel wall, estimated the cryogen boiloff, 
and showed the effects the ullage condition has on 
pressure-vessel temperatures. The degree of ullage 
mixing, location of the applied high-temperature pro- 
file, and the purge gas influence on insulation thermal 
conductivity had significant effects on the thermal be- 
havior of the Generic Research Cryogenic Tank. In ad- 
dition to analysis results, a description of the Generic 
Research Cryogenic Tank and the role it will play in 
future thermal structures and transatmospheric vehicle 
research at the NASA Dryden Flight Research Facility 
is presented. 
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Nomenclature 

Cp 

specific heat, BTU/lbm °R 

Gr 

Grashof number = g/3( Tw — Tv) x 3 / v 2 

gh 2 

gaseous hydrogen 

GRCT 

Generic Research Cryogenic Tank 

9 

gravitational constant, 32.2 ft/sec 2 

HA85M 

liquid Hydrogen test case, high temper- 
ature applied uniformly to All heat 
shield quadrants, 85-percent fill 
level, Mixed ullage condition 

HA85S 

liquid Hydrogen test case, high temper- 
ature applied uniformly to All heat 
shield quadrants, 85-percent fill 
level. Stratified ullage condition 

HB85M 

liquid Hydrogen test case, high temper- 
ature applied to the Bottom heat 
shield quadrant, 85-percent fill level, 
Mixed ullage condition 

HB85S 

liquid Hydrogen test case, high temper- 
ature applied to the Bottom heat 
shield quadrant, 85-percent fill level, 
Stratified ullage condition 

HT85M 

liquid Hydrogen test case, high temper- 
ature applied to the Top heat shield 
quadrant, 85-percent fill level. Mixed 
ullage condition 

HT85S 

liquid Hydrogen test case, high temper- 
ature applied to the Top heat shield 
quadrant, 85-percent fill level. 
Stratified ullage condition 

h 

convective heat transfer coefficient, 
BTU/ft 2 sec °R 



k thermal conductivity, BTU/ft sec °R 

LHSTF Liquid Hydrogen Structural Test Facility, 
Dryden Flight Research Facility, 
Edwards, CA 

LH 2 liquid hydrogen 

NIST National Institute of Standards and 
Technology, Boulder, CO 

Pr Prandtl number = y.c p /k 

Ra Rayleigh number = Gr Pr 

Re Reynolds number = ux/v 

S circumferential surface distance, ft 

T temperature, °R 

TAV transatmospheric vehicle 

u characteristic velocity, ft/sec 

x characteristic length, ft 

Y vertical distance, ft 

ft volumetric thermal expansion coeffi- 

cient, 1/°R 

/x dynamic viscosity, Ibm sec/ft 2 

v kinematic viscosity, ft 2 /sec 

1- D one-dimensional 

2- D two-dimensional 

Subscripts 

V vapor 

W wall 

Introduction 

Transatmospheric vehicles (TAVs) such as the Na- 
tional Aerospace Plane will require a fuselage which 
can withstand high aerodynamic heating while provid- 
ing an insulation system. This insulation system must 
reduce the heat load imposed on liquid hydrogen con- 
tained within onboard fuel tanks. Material degrada- 
tion, which occurs at the elevated surface temperatures 
associated with aerodynamic heating, restricts or dis- 
qualifies the use of many standard cryogenic insulat- 
ing materials, such as closed-cell foams or vacuum- 
jacketed multi-layer insulations. Thermal gradients 
which develop within the walls of the fuel tank can 
lead to high thermal stresses that affect tank integrity. 
Therefore, the development of new insulating systems 
for cryogenic fuel tanks and the validation of tank in- 


tegrity over a wide range of flight conditions will re- 
quire extensive testing. 

Tankage systems for TAVs have a significant impact 
on the overall vehicle design and have been the subject 
of several experimental test programs/ 1-3 ) These test 
programs helped to design, fabricate, and obtain ex- 
perimental validation of liquid hydrogen tankage ap- 
plicable to vehicles in hypersonic environments. Be- 
cause of the complex thermal interactions between the 
cryogenic fuel and the tank structure, the numerical 
simulation and optimization analysis of tank designs 
have also been an integral part of experimental test 
programs/ 4 

Personnel at the NASA Dryden Flight Research Fa- 
cility in Edwards, California are currently involved in 
the design of the Liquid Hydrogen Structural Test Fa- 
cility (LHSTF), to be completed in late 1993. When 
completed, the LHSTF will be able to test various full- 
scale and sub-scale flight vehicle components in simul- 
taneous cryogenic and high-temperature environments 
combined with mechanical loads. The LHSTF design 
consists of a large test cell for evaluating the perfor- 
mance and integrity of proposed TAV fuel tanks and 
associated insulation systems. In addition, the LHSTF 
site layout will provide for future capabilities includ- 
ing an actively cooled panel and turbomachinery test 
capability, altitude simulation, and full-scale vehicle 
fuselage and integrated systems tests. 

In preparation for cryogenic test operations, person- 
nel from NASA Dryden and PRC Inc. (formerly Plan- 
ning Research Corporation) have designed the Generic 
Research Cryogenic Tank (GRCT) as the first test arti- 
cle scheduled for testing in the LHSTF. As a research 
tank, the GRCT was designed to qualitatively simulate 
the thermal response of a TAV fuel tank exposed to the 
environment of hypersonic flight. The GRCT was de- 
signed to be a sturdy test article capable of withstand- 
ing a variety of operational and research tests. 

In studying the GRCT, NASA personnel will gain 
experience in operating, testing, and analyzing struc- 
tures in simultaneous cryogenic and high-temperature 
environments. Operationally, the GRCT will allow 
NASA personnel to learn how to handle cryogens 
and the associated equipment required for conducting 
cryogenic tests. Test operations with the GRCT will 
help develop and verify instrumentation capable of op- 
erating in both cryogenic and high-temperature envi- 
ronments, and help evaluate the thermal performance 
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of various insulation systems. To augment analysis 
efforts, the GRCT will provide test data for refining 
numerical models developed for simulating the ther- 
mal response of cryogenic tanks. Initial tests of the 
GRCT containing liquid nitrogen will be conducted 
in the high bay of the NASA Dryden Thermostruc- 
tures Research Facility with subsequent liquid hydro- 
gen tests conducted in the LHSTF. In addition, the 
GRCT will be used to perform the LHSTF integrated 
systems test before testing large and expensive TAV 
cryogenic fuel tanks. 

Numerical models were required to evaluate the 
thermal performance of the GRCT and to answer vari- 
ous design questions about insulation thickness, purge 
gas effects, and temperature gradients within the wall 
of the pressure vessel. To assist in the design and 
evaluation of the GRCT thermal performance, one- 


dimensional (1-D) and two-dimensional (2-D) ther- 
mal models were created using the Systems Improved 
Numerical Differencing Analyzer and Fluid Integra- 
tor (SINDA’85/FLUINT). (6) The SINDA’85/FLUINT 
uses a finite-difference solution method for analyzing 
thermal-fluid systems. 

This paper describes the components of the GRCT, 
the development of the 1-D and 2-D thermal models, 
and the response of the numerical models to several 
GRCT liquid hydrogen test scenarios. Transient re- 
sults for the 1-D and 2-D thermal models are presented 
with a discussion of the effects the numerical results 
had on the design of the GRCT. 

Overview of the GRCT 

Figure 1 shows a perspective view of the GRCT, 
suspended below a steel support structure, without 
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the piping and heat lamps required for testing. Fig- sel. This purge region serves two purposes. First, 

ure 2 shows a cut-away view of the GRCT along the since some anticipated liquid hydrogen test scenarios 



Section A-A 


Fig. 2 Cut-away view of the GRCT and a section view of the test section (not to scale). 


centerline and a section view through the 10-fit center 
test section. As shown in Figs. 1 and 2, the GRCT 
consists of a cylindrical stainless-steel pressure ves- 
sel (5-ft diameter by 15-fit length by 5/16-in. thick) 
surrounded by 3 in. of fibrous alumina-silica ceramic 
insulation blankets (8 lbm/ft 3 density) surrounded by 
a thin Inconel® heat shield (0.030-in. thick). A 
purge jacket of 0.005-in. thick nickel foil is located 
within the insulation at 1.5 in. from the pressure ves- 
sel. Purge gas, helium for liquid hydrogen tests and 
nitrogen for liquid nitrogen tests, will be pumped into 
the end bell sections of the GRCT and channeled into 
the inner 1.5 in. of insulation next to the pressure ves- 

® Inconel is a registered trademark of Huntington Alloy Products 
Division, International Nickel Company, Huntington, WV. 


will be conducted in an air atmosphere within the 
LHSTF, a helium purge region reduces the possibil- 
ity of a hydrogen leak forming a combustible mixture 
with air. Second, the purge region eliminates the liq- 
uefaction of air or nitrogen (the two LHSTF test cell 
atmospheres) within the insulation which degrades in- 
sulation performance and creates potential safety and 
maintenance problems. 

Allowing for an adequate ullage (the unfilled por- 
tion of a container) for pressure relief, the maximum 
fill level for the GRCT will be from 85 to 90 percent 
of the total pressure vessel volume of 267 ft 3 . A fill- 
drain line to the pressure vessel simulates the liquid 
cryogen outflows required during TAV flight profiles 
for engine and cooling requirements. The fill-drain 
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line has been sized to provide a maximum cryogen out- 
flow rate of 2 lbm/sec. A pressurization line was incor- 
porated in the pressure vessel design to provide tank 
pressure maintenance during cryogen outflow. Vent 
lines, with back-pressure regulators, are used to con- 
trol the pressure within the pressure vessel and allow 
a maximum pressure of 45 psia. Instrumentation ac- 
cess is provided for temperature, pressure, and liquid 
level measurements within the pressure vessel. As a 
design requirement, the fill-drain, pressurization, vent, 
and instrumentation penetrations are confined to the 
hemispherical ends of the GRCT to provide a uniform 
cylindrical test section clear of penetrations. 

During test operations, clamshell quartz lamp 
heaters will be placed around the suspended GRCT. 
The quartz lamps will heat the outer heat shields and 
provide the temperature load on the GRCT. Figure 3 
shows the proposed heating profiles to be applied to the 


fects on the ullage, the high-temperature profile will 
be applied to the GRCT upper heat shield quadrant 
(hot-top) while the low-temperature profile (peak tem- 
perature of 1260 °R) is applied to the lower quadrant. 
To examine high-temperature effects on the liquid re- 
gion, the profiles will be reversed (hot-bottom). Dur- 
ing nonuniform heating, the side quadrants of the heat 
shield will follow a heating profile composed of an av- 
erage of the high- and low-temperature profiles. 

Description of the Thermal Models 

The GRCT design provides a test capability with 
either liquid nitrogen or liquid hydrogen. The ther- 
mal models developed to simulate the response of the 
GRCT also simulate either liquid nitrogen or hydro- 
gen. However, analyzing the GRCT design focused on 
the thermal response associated with liquid hydrogen 



Fig. 3 Proposed heating profiles applied to the heat shields of the GRCT. 


heat shields of the GRCT. These temperature pro- 
files represent the temperature loads associated with 
TAV flight trajectories. A baseline even-heating test 
case will consist of the high-temperature profile (peak 
temperature of 1960 °R) uniformly applied to the heat 
shield. Tb investigate the effects of nonuniform heat- 
ing associated with TAV ascent and descent flight tra- 
jectories, the two heating profiles shown in Fig. 3 will 
be applied nonuniformly. To examine temperature ef- 


testing and the results presented are confined to liquid 
hydrogen test conditions. 

One-Dimensional Thermal Model 

The GRCT 1-D thermal model was developed to de- 
termine the required insulation thickness and the lo- 
cation of the purge liner. The 1-D model represented 
a radial section of the GRCT from the liquid hydro- 
gen to the heat shield and accounted for conduction. 
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convection, and radiation heat transfer. Figure 4 is a 
schematic of the node and conductor layout used for 


ture of hydrogen for 45 psia. The amount of boiloff 
produced during each time step was calculated by 
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Fig. 4 Representative schematic of the node and conductor layout for the 1-D thermal model. 


the 1-D thermal model. A total of 36 nodes, each 
with a constant 1 ft 2 area, were used in the model. 
The thermal conductivity of the Inconel heat shield, 
insulation, and stainless-steel pressure vessel was tab- 
ulated as functions of temperature/ 7 - 8 ) Heat transfer 
through the insulation was modeled entirely by con- 
duction, while conduction and radiation heat transfer 
were included between the heat shield and the outer in- 
sulation layer. A constant effective emissivity of 0.53 
and a view factor of 1 .0 were assumed for the radiation 
heat transfer between the heat shield and the insulation. 

All computational test scenarios were examined as- 
suming a constant tank pressure of 45 psia. Liq- 
uid hydrogen was assumed to be in continuous con- 
tact with the pressure-vessel wall with nucleate boil- 
ing defined as the mode of heat transfer for the 
wall-to-liquid interface. The heat transfer coeffi- 
cient between the pressure- vessel wall and the liq- 
uid hydrogen was estimated from the correspond- 
ing Kutateladze correlation for nucleate pool boil- 
ing and tabulated as a function of the temperature 
difference between the pressure vessel and the liq- 
uid cryogen/ 9 ) Nucleate boiling was assumed to pro- 
duce a well mixed liquid region, which allowed the 
liquid node to be held at the saturation tempera- 


dividing the heat transferred into the liquid node by the 
enthalpy of vaporization for hydrogen (193 BTU/lbm). 

A problem with the insulation thermal conductivity 
was identified during the 1-D thermal modeling work. 
Insulating materials are often not well characterized 
for applications involving large temperature gradients. 
High-temperature conductivity values are often mea- 
sured using a small imposed temperature gradient, but 
materials used within the GRCT will be exposed to 
gradients in excess of 600 °R per inch. Under these 
conditions, fiber-to-fiber radiation can contribute to the 
overall apparent thermal conductivity of the insulation. 
The impact of purge gas (nitrogen or helium) on the ap- 
parent thermal conductivity must also be considered. 
Gas conduction is the dominant heat transfer mech- 
anism for fibrous insulation systems at atmospheric 
pressure/ 10 ) Therefore, because helium is more con- 
ductive than air, the characteristics of the apparent 
thermal conductivity for the selected fibrous insulation 
will change in the presence of a helium purge. A pro- 
gram is currently underway with the National Institute 
of Standards and Technology (NIST) to measure the 
apparent conductivity of the GRCT fibrous alumina- 
silica insulation at high temperature gradients coupled 
with helium and nitrogen purge gases. The addition 
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of the resulting apparent conductivity data from NIST 
will enhance the prediction capability of the thermal 
models developed for the GRCT. 

The NIST data were not available for this study, 
hence, the manufacturer’s thermal conductivity data 
for the insulation was used within the GRCT thermal 
models. However, since the insulation conductivity 
data were obtained by the manufacturer in an air atmo- 
sphere, the effects of a helium purge gas on the insula- 
tion thermal conductivity had to be modeled. To esti- 
mate the helium purge effects on the fibrous insulation, 
tabulated data for the thermal conductivity of perlite, 
a silica powder, for densities from 6 to 9 lbm/ft 3 and 
for various interstitial gases at one atmosphere were 
examined/ 11 ) These data showed that the ratio of the 
perlite thermal conductivity containing helium to that 
containing nitrogen was approximately three for all 
densities. Therefore, to simulate the effects of a he- 
lium purge gas within the GRCT insulation, the man- 
ufacturer’s thermal conductivity data was multiplied 
by three within the purge region (inner 1.5 in.) and re- 
mained unchanged in the outer insulation region (outer 
1.5 in.). 

The Two-Dimensional Thermal Model 

The 2-D thermal model of the GRCT was created 
to examine the temperature gradients developed within 
the pressure-vessel wall, refine the calculation of cryo- 
gen boiloff, and characterize the thermal behavior of 
the ullage. The 2-D model represented a 1 -ft wide 
portion of the GRCT cylindrical test section and mod- 
eled the hydrogen liquid and ullage regions. Figure 5 
is a schematic of the node and conductor layout used 
for the 2-D thermal model composed of 413 nodes. 
Twenty-six circumferential locations were defined be- 
tween the pressure vessel and the heat shield bound- 
aries. At each circumferential location, 15 nodes were 
used to model the radial and circumferential heat trans- 
fer from the pressure-vessel wall, through the insu- 
lation, and to the heat shield. Each node within this 
region had a cross-sectional area and volume propor- 
tional to its radial location. The hydrogen contained 
in the GRCT was modeled by dividing the pressure- 
vessel cross-sectional area into 12 horizontal sections 
from top to bottom, with the uppermost horizontal sec- 
tion containing 3 vapor nodes. When the pressure ves- 
sel was 85-percent full, the top 5 fluid nodes were va- 
por and the remaining nodes were defined as liquid (a 
total of 18 nodes in liquid). Based on results of the 


1-D model, only one radial node was defined within 
the pressure- vessel wall at each circumferential loca- 
tion since the radial temperature gradients were small. 
The 2-D thermal model assumed a constant tank op- 
erating pressure of 45 psia and assumed the liquid hy- 
drogen was well mixed and remained at the prescribed 
saturated liquid temperature. 

The 2-D model simulated the effects of hydrogen 
boiloff on the ullage, in addition to conduction, con- 
vection, and radiation heat transfer. The specific heat 
of the hydrogen vapor nodes and the boiloff gas were 
tabulated as a function of temperature/ 12) As with the 
1-D model, the wall-to-liquid heat transfer coefficient 
for the liquid hydrogen cases was based on the nucle- 
ate pool boiling correlation of Kutateladze however, 
for stability reasons, the correlation was tabulated as a 
function of time rather than temperature difference. 

Both free and forced convection correlations were 
considered in modeling the wall-to-vapor heat transfer 
within the ullage, however, an estimate was made to 
determine if either was a dominant convection mode 
or if mixed convection best described the ullage heat 
transfer. The flow velocity developed by liquid hydro- 
gen boiloff produced laminar-flow forced convection, 
while the temperature difference between the vapor 
and the pressure-vessel wall produced free convection. 
An aid to determine the dominant convective mode 
of heat transfer was the ratio of the Grashof number 
to the Reynolds number squared (Gr/Re 2 ) which pro- 
vides a measure of the ratio of buoyancy forces (free 
convection) to inertial forces (forced convection)/ 13) 
A flow is considered dominated by free convection 
if Gr/Re 2 ;» 1 and dominated by forced convection 
if Gr/Re 2 < 1. The liquid hydrogen 85-percent fill 
level test case yielded a Gr/Re 2 ratio on the order of 
1 x 10 5 which suggested that the wall-to-vapor heat 
transfer was dominated by free convection. Deter- 
mining whether the free convective flow was laminar 
or turbulent was deduced from the Rayleigh number 
(Ra), which is a measure of the ratio of the buoy- 
ancy forces to the viscous forces. A free convec- 
tive flow is considered laminar if Ra < 10 9 and con- 
sidered turbulent if Ra > 10 9 . The Ra for the liq- 
uid hydrogen 85-percent fill level test case was on 
the order of 10 12 to 10 14 . As a result, the wall-to- 
vapor heat transfer within the ullage was modeled as 
turbulent-free convection and calculated from the tem- 
perature difference between the pressure-vessel wall 
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Fig. 5 Representative schematic of the node and conductor layout for the 2-D thermal model. 


and the hydrogen vapor. The form of the turbulent- 
free convection correlation^ 4 ^ used in the 2-D model 
is given by 

h = 0.13-(GrPr) 1/3 

x 

where k and Pr were the thermal conductivity and 
Prandtl number of hydrogen, and x was a character- 
istic length (the pressure vessel diameter of 5 ft). 

Test cases representing a stratified or well mixed ul- 
lage condition were examined by modifying the vapor- 
to-vapor heat transfer coefficient within the ullage. To 
verify the GRCT design under worst case conditions, 


the ullage vapor-to-vapor heat transfer coefficient was 
modeled by gaseous conduction (corresponding to heat 
transfer coefficients of 4 x 10~ 6 to 2 x 10 -5 BTU/ft 2 
hr °R) which produced a high degree of ullage strati- 
fication. To simulate a well mixed ullage, the vapor- 
to-vapor heat transfer coefficients (based on gaseous 
conduction) were multiplied by a factor of 1 x 10 6 
to yield heat transfer coefficients on the same order of 
magnitude as the wall-to- vapor free convection coeffi- 
cients (approximately 9 to 20 BTU/ft 2 hr °R). 

Table 1 shows the nomenclature used to identify the 
conditions of the various liquid hydrogen test cases ex- 
amined using the 2-D thermal model. 
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Table 1. Computational test matrix examined with the 2-D thermal model of the GRCT. 


Arrangement of the heating profiles 
applied to the heat shield quadrants 

Test case identification 
(Liquid hydrogen 85-percent fill level) 
Stratified ullage Well mixed ullage 

Nonuniform heating 
High-temperature profile-top quadrant 
Low-temperature profile-bottom quadrant 

HT85S 

HT85M 

Uniform heating 

High-temperature profile-all quadrants 

HA85S 

HA85M 

Nonuniform heating 

High-temperature profile-bottom quadrant 
Low-temperature profile-top quadrant 

HB85S 

HB85M 


Each test case identified in Table 1 describes the 
type of cryogen modeled, the location of the high- 
temperature profile, the fill level, and the condition 
of the ullage. For example, the HT85S identification 
translates to liquid hydrogen within the pressure ves- 
sel, the high-temperature profile applied to the top heat 
shield quadrant, an 85-percent fill level, and a stratified 
ullage. All of the other test cases in Table 1 follow a 
similar nomenclature format. 

Computational Results 
One-Dimensional Results 

The 1-D thermal model represented the heat trans- 
fer associated with the pressure vessel in contact with 
liquid hydrogen and did not simulate the heat transfer 
within the ullage of the GRCT. Therefore, all boiloff 
calculations within the 1-D results section represent an 
estimate of the total pressure-vessel boiloff rate based 
on a 100-percent fill level. The boiloff rate for a 100- 
percent fill level was calculated by multiplying the 1-D 
model boiloff rate by the total pressure- vessel internal 
surface area of 236 fit 2 . 

Insulation Thickness 

The effect of insulation thickness on the thermal 
response of the GRCT was examined using the 1-D 
model. The design criteria for the GRCT required the 
insulation system to provide a steady-state pressure- 
vessel wall heat flux of approximately 30 BTU/ft 2 hr, 
simulating the ground hold condition of a possible 
TAV design. In addition, design criteria for the peak 
wall heat flux required at least an order of magni- 
tude increase from steady-state conditions (to approx- 
imately 300 BTU/ft 2 hr) within the 3000-sec applied 
heating period. 


Using the manufacturer’s thermal conductivity data, 
Fig. 6 shows the pressure-vessel wall heat flux pre- 
dicted for liquid hydrogen and insulation thicknesses 
of 2 to 6 in. The wall heat flux associated with 
3 in. of insulation produced a steady-state wall flux 
of 22 BTU/ft 2 hr, corresponding to a steady-state 
boiloff of 27 lbm/hr, and yielded a peak heat flux of 
153 BTU/ft 2 hr which equates to a peak boiloff rate of 
187 lbm/hr. The 2-in. insulation thickness produced 
a steady-state wall heat flux of 33 BTU/ft 2 hr, how- 
ever, it produced an excessive peak heat flux and corre- 
sponding boiloff rate (384 BTU/ft 2 hr and 469 lbm/hr) 
within the 3000-sec heating period. The 4- to 6-in. 
range of insulation thicknesses yielded low steady- 
state wall heat fluxes (16 to 11 BTU/ft 2 hr) and pro- 
duced virtually no thermal response within the tran- 
sient heating period. Because of the anticipated in- 
crease in wall heat flux associated with using a helium 
purge, the 2-in. insulation option was eliminated and 
3 in. of insulation was selected for the region between 
the pressure vessel and the heat shield. 

Location of the Purge Liner 

At one atmosphere, air liquefies at 142 °R and 
any air within the insulation below this tempera- 
ture liquefies, degrading insulation performance and 
creating potential safety and maintenance problems. 
To eliminate the air liquefaction problem, the nickel 
foil purge jacket was placed within the insulation 
at 1.5 in. to ensure the insulation temperature out- 
side the jacket remained above 142 °R. The steady- 
state hold condition for the GRCT (that is, the 
GRCT filled with liquid cryogen awaiting the start 
of a test) produced the lowest insulation tempera- 
tures corresponding to the highest probability of 


9 




Fig. 6 Pressure vessel wall heat flux as a function of time from the 1-D thermal model with liquid hydrogen and 
insulation thicknesses from 2 to 6 in. 


liquefying air. The liquid hydrogen filled GRCT pro- 
duced the largest temperature difference between the 
liquid cryogen and the ambient environment, thereby 
defining the radial location of the purge jacket. 

Figure 7 shows the steady-state and maximum tran- 
sient temperature distributions through the fibrous 


insulation with the pressure vessel containing liquid 
hydrogen. The helium purge was confined to the inner 
1.5 in. of insulation next to the pressure vessel. The 
maximum transient curve shows the peak tempera- 
tures that occur within the insulation, while the steady- 
state curve defines the minimum insulation tempera- 
tures associated with the ground hold condition. This 



Fig. 7 Steady-state and maximum transient insulation temperature distribution as a function of radial location for 
a 3-in insulation blanket with a helium purge. 
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figure shows that with the purge liner located 1.5 in. 
from the pressure vessel, the temperature of the un- 
purged region will remain above the air liquefaction 
temperature. 

Helium Purge Gas Effects 

Figure 8 is a comparison of the pressure-vessel wall 
heat flux for 3 in. of insulation with and without a 
helium purge. The addition of the helium purge gas 



Fig. 8 The 1-D thermal model comparison of pres- 
sure vessel wall heat fluxes for a 3-in insulation blan- 
ket with and without helium purge as a function 
of time. 

increased the steady-state wall heat flux of the GRCT 
from 22 BTU/ft 2 hr to 33 BTU/ft 2 hr, and increased 
the steady-state boiloff from 27 lbm/hr to 40 lbm/hr. 
The peak transient wall heat flux increased from 
153 BTU/ft 2 hr to 312 BTU/ft 2 hr corresponding to 
an increase in peak transient boiloff from 187 lbm/hr 
to 381 lbm/hr. The helium purge gas not only affected 
the quantitative values of the wall heat flux and boiloff 
but also the transient behavior of the system. The he- 
lium purge gas decreased the thermal resistance of the 
inner 1.5 in. of insulation which caused the pressure 
vessel to respond more quickly to the imposed thermal 
environment. 

A flexible GRCT design was needed because of 
uncertainty in the helium purge effects on the insu- 


lation conductivity. Two methods of controlling the 
pressure-vessel wall heat flux values have been iden- 
tified in case the helium purge effects produce higher 
pressure-vessel wall heat flux values than calculated by 
the thermal model. First, to reduce the wall heat flux 
values, the heat shield has been designed to accommo- 
date an additional inch of insulation, yielding a total 
of 4 in. of insulation if required. Second, the heating 
profiles may be altered to achieve a desired heat flux 
at the pressure- vessel wall. 

Two-Dimensional Results 

The 2-D thermal model simulated the heat transfer 
interaction between the ullage and liquid regions of the 
pressure vessel. Therefore, the total pressure-vessel 
boiloff rate presented in the 2-D results section for an 
85-percent fill level was calculated by multiplying the 
2-D model boiloff rate by the ratio of the total pressure- 
vessel inner surface area (236 ft 2 ) to the 2-D model 
pressure-vessel surface area (15.7 ft 2 ). In comparing 
1 -D and 2-D boiloff rates, the 1 -D boiloff rate was mul- 
tiplied by the total wetted surface area of the pressure 
vessel for an 85-percent fill level (174 ft 2 ). 

Liquid Hydrogen Boiloff 

Figure 9 compares the liquid hydrogen boiloff for 
the 1-D and 2-D thermal models and the HA85S 



Fig. 9 Comparison of liquid hydrogen boiloff as a 
function of time for the 1-D and 2-D thermal models 
and the HA85S test conditions. 
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test conditions. Properly scaled for an 85-percent fill 
level, the 1-D thermal model yielded a transient boiloff 
rate similar to the 2-D thermal model. For the HA85S 
test case, the 2-D thermal model calculated a steady- 
state boiloff of 32 lbm/hr and a peak boiloff of 305 
lbm/hr, compared with 30 lbm/hr and 281 lbm/hr for 
the 1-D model. Since the additional heat transfer be- 
tween the ullage and liquid regions was included, the 
2-D thermal model provided a more realistic estimate 
of the transient hydrogen boiloff. 

Figure 10 shows a comparison of the liquid hy- 
drogen boiloff as a function of time for the HT85S, 
HA85S, and HB85S test cases. This figure demon- 


duced boiloff. With the high-temperature profile in the 
top quadrant, a lower boiloff was produced because 
the hydrogen vapor simply absorbed the heat and in- 
creased in temperature. 

Pressure- Vessel Wall Temperatures 

The maximum wall temperature for the 6 test cases 
at the 85-percent fill level was determined by ex- 
amining the transient behavior of the pressure-vessel 
wall temperature in the ullage. Figure 11 shows the 
pressure-vessel wall temperature as a function of time 
for the HT85S test case and for computational nodes 
100, 400, and 500. The peak wall temperature of 



Fig. 10 Comparison of liquid hydrogen boiloff as a function of time for different applied temperature profile 
arrangements with the 2-D thermal model. 


strates the effect the high- temperature profile loca- 
tion had on the cryogen boiloff rate. The high- 
est liquid hydrogen boiloff rate occurred when the 
high-temperature profile was applied to all heat shield 
quadrants which yielded the maximum wetted sur- 
face area in contact with high wall heat fluxes. For 
the nonuniformly applied heating profiles, HB85S 
yielded a higher boiloff than HT85S. With the high- 
temperature profile on the bottom quadrant, the liquid 
hydrogen readily absorbed the wall heat flux and pro- 


198 °R occurred at the highest point in the pressure 
vessel (node 100) at 3600 sec. However, the peak 
temperatures for each node within the pressure-vessel 
wall did not all occur at the same time. For exam- 
ple, the peak temperature of node 400, 75 °R, occurred 
at 2800 sec. This time difference for the peak wall 
temperatures was caused by the geometric arrange- 
ment of the vapor nodes and the resulting boiloff mass 
distribution. 
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Fig. 11 Pressure vessel wall temperatures as a function of time for the HT85S test case and computational nodes 
100, 400, and 500 with the 2-D thermal model. 


The HT85S test case (hot-top profile) yielded the 
peak pressure-vessel wall temperature because the 
high-temperature profile on the top heat shield quad- 
rant was combined with the lowest boiloff rate (see 
Fig. 10). Figure 12 shows the transient temperature 
of node 100 for the 3 locations of the applied high- 


temperature profile. As shown in Fig. 12, concentrat- 
ing the high-temperature profile on the top heat shield 
quadrant yields the lowest boiloff rate, resulting in 
less cooling available for the wall and ullage regions, 
which yielded the highest wall temperatures. Con- 
versely, concentrating the high-temperature profile on 



Fig. 12 Pressure vessel wall temperature for node 100 as a function of time for different applied temperature 
profile arrangements with the 2-D thermal model. 
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the bottom heat shield quadrant yielded a moderate 
boilofT rate, which coupled with the low-temperature 
profile in the ullage, produced lower overall wall tem- 
peratures. 

The liquid hydrogen acted as a thermal sink which 
maintained the wetted pressure-vessel wall tempera- 
tures slightly above the saturation temperature of the 
hydrogen for all of the applied temperature profile ar- 
rangements examined. The magnitude of the wall- 
to-liquid heat transfer coefficient for liquid hydrogen 
did not substantially affect the value of the wetted 
pressure-vessel wall temperatures. 

While examining the ullage heat transfer, various 
constant values for the wall-to-vapor convective heat 
transfer coefficient (from 0.5 to 100 BTU/ft 2 hr °R) 
were considered. This sensitivity analysis showed 
that the wall-to-vapor heat transfer coefficient had 
a minor effect on the pressure-vessel wall and fluid 
temperatures. Conversely, the magnitude of the 
vapor-to-vapor heat transfer coefficient (affecting the 
degree of mixing) greatly influenced the resulting 
pressure- vessel wall and fluid temperatures. Of 
all the modes of heat transfer occurring within the 


pressure vessel, the condition of the ullage va- 
por (well mixed or stratified) had the most pro- 
nounced effect on the pressure-vessel wall temper- 
atures. A stratified ullage yielded higher pressure- 
vessel wall temperatures and circumferential tempera- 
ture gradients while a well mixed ullage yielded lower 
wall temperatures and circumferential temperature 
gradients. 

Figure 13 shows the pressure- vessel wall tempera- 
ture distribution for the HT85S and the HT85M test 
cases. For the HT85S test case, the pressure-vessel 
wall temperatures are shown at 0 sec (steady-state) 
and at 3600 sec (peak temperature). The pressure- 
vessel wall temperatures decreased circumferentially 
from S = 0 ft (node 100, top of the pressure ves- 
sel) to S = 2.6 ft (node 500) with the heat flow- 
ing from the hot upper wall region to the cold wall 
region in contact with the liquid hydrogen. Figure 
13 also provides an estimate of the circumferential 
temperature gradients that exist within the pressure- 
vessel wall. Comparing the temperature distributions 
for the HT85S and HB85M test cases shows the de- 
crease in the pressure-vessel wall temperatures and 
temperature gradients when the ullage was well mixed. 



Fig. 13 Steady-state and peak pressure vessel wall temperatures as a function of circumferential location for 
stratified and mixed ullage conditions with the 2-D thermal model. 
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Vapor Temperatures 

Figure 14 shows the vertical temperature distri- 
bution of the fluid and vapor nodes at selected 
times for the HT85S and HT85M test cases. 


function of the location of the applied heating pro- 
files. The heating profiles influenced the amount of 
cryogen boiloff which then influenced the vapor tem- 
peratures and the corresponding wall temperatures. 



Vertical position, Y, ft 910129 

Fig. 14 Liquid and vapor temperatures as a function of vertical position for a stratified and well mixed ullage 
condition with the 2-D thermal model. 


Node 3000 (the top vapor node) reached a peak vapor 
temperature of 188 °R at 3600 sec, which was 10 °R 
cooler than the adjacent wall nodes (see Fig. 11). The 
temperature gradient associated with the stratified ul- 
lage and its change over time is readily apparent in 
Fig. 14. For a well mixed ullage condition there was 
virtually no temperature gradient through the vapor re- 
gion as shown in Fig. 14. The time dependence of the 
fluid and vapor node temperatures was qualitatively 
similar to the time dependence of the wall tempera- 
tures shown in Figs. 11 and 12. 

Pressure Vessel Circumferential Temperature 
Gradients 

The circumferential temperature gradients devel- 
oped within the pressure-vessel wall were a strong 


The 2-D thermal model was used to examine the peak 
thermal gradients developed within the pressure-vessel 
wall for the 85-percent fill level and determine if they 
were below the design criteria of 112 °R/ft prescribed 
by the structural analysis done by PRC personnel. 

Table 2 shows the peak and average circumferential 
temperature gradients developed within the pressure- 
vessel wall for the different applied heating profiles 
with stratified and well mixed ullages. The circum- 
ferentially averaged temperature gradients were calcu- 
lated from the pressure-vessel temperature gradients in 
contact with the ullage at the time the peak gradients 
occurred (typically 2500 to 5000 sec). 
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Table 2. Peak and average circumferential temperature 
gradients within the pressure-vessel wall for stratified 
and well mixed ullages. 


Test cases 

Pressure vessel circumferential 
temperature gradients (°R/ft) 


Peak 

Average 

Stratified ullage 

HT85S 

116.1 

58.2 

HA85S 

99.9 

44.8 

HB85S 

54.3 

25.4 

Well mixed ullage 

HT85M 

21.0 

6.6 

HA85M 

25.4 

6.4 

HB85M 

13.0 

3.2 


The peak gradient always occurred as a spike that 
was significantly larger than the rest of the gradients 
in the ullage. These temperature gradient spikes were 
the result of several simplifying assumptions built 
into the 2-D model, including the degree of vapor-to- 
vapor coupling and the distribution of boiloff vapor 
(for cooling) available to each tank wall node. Cir- 
cumferentially averaging the ullage temperature gra- 
dients reduced the effects of the simplifying model 
assumptions. Therefore, the averaged values were 
considered to be more appropriate for design assess- 
ment. The average temperature gradients for the 
stratified ullage case were considerably less than the 


design criteria. A well mixed ullage yielded average 
temperature gradients which were considerably lower 
than for a stratified ullage and were well below the de- 
sign criteria. 

Pressure- Vessel Wall Heat Flux 

Table 3 shows the peak and average pressure-vessel 
wall heat fluxes with liquid hydrogen for stratified and 
well mixed ullages. The average wall heat fluxes were 
calculated from the nodes in contact with liquid or va- 
por at the time the corresponding peak wall heat fluxes 
occurred. 


Table 3. Peak and average pressure- vessel wall heat fluxes with liquid hydrogen 
for stratified and well mixed ullages. 


Test cases 

Wall-to-vapor heat fluxes 
(BTU/ft 2 hr) 

Wall-to-liquid heat fluxes 
(BTU/ft 2 hr) 


Peak 

Average 

Peak 

Average 

Stratified ullage 

HT85S 

348.9 

241.9 

208.0 

164.5 

HA85S 

378.2 

309.3 

351.4 

342.6 

HB85S 

140.7 

112.5 

341.3 

244.0 

Well mixed ullage 

HT85M 

340.6 

318.5 

200.2 

163.5 

HA85M 

342.2 

340.7 

350.7 

347.0 

HB85M 

138.6 

114.3 

343.9 

236.6 
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The location and magnitude of the peak wall heat 
flux within the ullage depended on several factors, in- 
cluding the applied high-temperature profile location 
and the geometric arrangement of the vapor nodes af- 
fecting boiloff mass distribution. Averaging the ullage 
temperature gradients reduced the effects of the simpli- 
fying model assumptions. Consequently, the average 
wall-to-vapor heat flux provided a more realistic indi- 
cation of the expected wall heat flux. For the hot-top 
and even heating test cases more heat was transferred 
into the ullage than for the hot-bottom test cases. Ul- 
lage mixing increased the wall-to-vapor heat flux for 
the hot-top and even heating cases, but had no effect 
on the hot-bottom cases. 

The magnitude of the peak wall-to-liquid heat flux 
was not dependent on the arrangement of the liquid 
nodes, but the location of the peak wall-to-liquid heat 
flux depended on the location of the applied high- 
temperature profile. For the hot-top and even heating 
cases, the peak heat flux occurred at the top wall node 
in contact with the liquid cryogen (node 500) because 
of the additional heat transferred from the ullage. For 
the hot-bottom heating cases, the peak wall-to-liquid 
heat fluxes occurred at the bottom of the pressure ves- 
sel. The condition of the ullage, whether stratified or 
well mixed, did not affect the wall-to-liquid heat fluxes 
within the pressure vessel. 

Concluding Remarks 

The one-dimensional and two-dimensional thermal 
models successfully analyzed the thermal behavior of 
the Generic Research Cryogenic Tank for several test 
cases. The insulation thickness around the pressure 
vessel was sized at 3 in. to provide a steady-state heat 
flux of 33 BTU/ft 2 hr and produced a peak transient 
heat flux of 312 BTU/ft 2 hr which occurred within the 
3000-sec heating period. The purge jacket was located 
1.5 in. within the insulation to eliminate the lique- 
faction of air which would degrade insulation perfor- 
mance. Large temperature gradients within the wall of 
the pressure vessel could lead to large thermal stresses, 
but the average circumferential temperature gradients 
of the test cases examined were well below the al- 
lowable design value of 112 °R/ft. Refining the cryo- 
gen boiloff calculation and subsequently the heat flux 
into the cryogen defined the thermal performance of 
the Generic Research Cryogenic Tank more accurately. 
Mixing in the ullage lowered vapor temperatures by as 


much as 140 °R and had a much greater effect on ther- 
mal behavior of the ullage than either heating profile 
location or wall-to-vapor heat transfer characteristics. 

The thermal analysis identified several characteris- 
tics that affected the thermal behavior of the Generic 
Research Cryogenic Tank. The characteristics of ul- 
lage mixing, the location of applied high-temperatures, 
and the helium-insulation apparent thermal conductiv- 
ity all influenced the behavior of the Generic Research 
Cryogenic Tank. Improvement of subsequent thermal 
simulations will be the focus of the thermal conductiv- 
ity tests at the National Institute of Science and Tech- 
nology and much of the future testing with the Generic 
Research Cryogenic Tank. 
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